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We study two microscopic models of topological insulators in contact with an s-wave supercon¬ 
ductor. In the first model the superconductor and the topological insulator are tunnel coupled via 
a layer of scalar and of randomly oriented spin impurities. Here, we require that spin-flip tunneling 
dominates over spin-conserving one. In the second model the tunnel coupling is realized by an array 
of single-level quantum dots with randomly oriented spins. It is shown that the tunnel region forms 
a TT-junction where the effective order parameter changes sign. Interestingly, due to the random 
spin orientation the effective descriptions of both models exhibit time-reversal symmetry. We then 
discuss how the proposed vr-junctions support topological superconductivity without magnetic fields 
and can be used to generate and manipulate Kramers pairs of Majorana fermions by gates. 

PACS numbers: 74.50.+r, 71.10.Pm, 74.45.+C 


Introduction. When two s-wave superconductors (SCs) 
are brought into contact via an insulator the energy of 
the system in equilibrium is minimized when the relative 
phase between the two superconducting order parameters 
vanishes. Interestingly, when the insulator is doped with 
magnetic impurities, it was shown by theory [T] and ex¬ 
periment [5] that spin-flip tunneling can induce an equi¬ 
librium ground state with a relative phase difference of tt 
between the superconducting order parameters, building 
up a so-called Josephson 7r-junction (JttJ). It was pre¬ 
dicted [5] and experimentally confirmed [1] that a JttJ 
can be generated by replacing the layer of magnetic im¬ 
purities by a ferromagnetic metal. A JttJ can also arise 
when two SCs are tunnel-coupled through an interme¬ 
diate resonant state in the presence of strong Coulomb 
interactions [5], as observed in a system of two SCs cou¬ 
pled by a quantum dot (QD) occupied by a single electron 
[5]. In recent experiments [7H5] it was demonstrated that 
superconductivity can also be proximity-induced in the 
helical edge states of a topological insulator (TI) mate¬ 
rial nm via coupling to an external s-wave SC. These 
experimental advances have also stimulated the theoreti¬ 
cal interest in Josephson junctions based on TIs [IMI]. 
Motivated by the existence of ordinary JttJs an impor¬ 
tant and immediate question is: Are there microscopic 
mechanisms allowing one to induce a superconducting 
order parameter in the helical edge states of the TI that 
is of opposite relative sign compared to the one of the ex¬ 
ternal s-wave SC, ideally without breaking time-reversal 
invariance (TRI)? In this work we answer this question 
in the affirmative. 

We propose two setups involving TIs in which such a 
TT-junction is shown to emerge. In the first setup the tun¬ 
nel coupling is realized by a thin insulating layer of scalar 
and magnetic impurities with randomly oriented spins. 
Here we require that spin-flip tunneling dominates over 
normal tunneling. In the second setup the tunnel cou¬ 
pling is realized by an array of single-level QDs, each of 
which is occupied by a single spin with random orienta¬ 


tion. Critically, the random orientation of spins preserves 
TRI in an effective description. We note that both setups 
can be realized by combining the already existing experi¬ 
ments on proximity-inducing superconductivity solely in 
the edge states of a TI [7HS] and the experiments on JttJs 
in SC-magnetic insulator-SC [5] and SC-QD-SC devices 
[B]. We note that the same setup could be assembled in 
the framework of strip of stripes models based on 

an array of coupled one-dimensional channels with spin- 
orbit interaction m- As a striking consequence we find 
that the proposed models for proximity-induced JttJs in a 
TI provide an alternative approach to engineer Kramers 
pairs of Majorana fermions (MFs) easily movable 

by gates. Remarkably, no magnetic fields are needed. 
More precisely we consider two TI samples that form a 
proximity-induced JttJs with respect to one another and 
allow for tunneling between them in the finite space re¬ 
gion, at the ends of which the MFs emerge. 

Josephson junetion models. In the first model we con¬ 
sider a bulk s-wave SC connected by a tunnel contact to 
the edge of a 2D TI, see Fig. Ba). The Hamiltonian of 
the system is given by 

Hi = Hbcs + Hti (1) 

- 1 -^ J dr dx [d'^(r) • Ti(r, a;)<i>(a;) -I-H.C.] , 

with the tunneling matrix Ti(r, a;) = Ti(r, a:)(l-|-r^)/2 — 
Ti(r,a:)(l - t"=)/ 2. Here, Hbcs = (l/2)/dr • 

[—{h‘^d^/2m + fx)T^ — AscCT^'r^]'I'(r) -|- H.c. is the BCS 
Hamiltonian of the SC, fi being the chemical poten¬ 
tial in the SC and m being the electron mass, and 
Hti = (1/2) / dx • {—ihvFO’^dx)^{x) + H.c.] is 

the Hamiltonian of the TI edge with the Fermi ve¬ 
locity vp- Without loss of generality, we assume 
that the superconducting order parameter A^c is posi¬ 
tive. The electron Nambu operator in the SC (TI) is 
given by d'(r) = (4'^(r), d';(r), d'|(r), «']](r)) [$(a:) = 
(R(a;),L(a;),Rl(x),Ll(a;))]. The Pauli matrices t“ ((t“) 
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with a = x,y,z act in particle-hole (spin) space The 
slowly-varying spin-up right (spin-down left) mover fields 
R(a;) [L(a;)] are defined around the Fermi points zLkp 
which in turn are determined by the position of the chem¬ 
ical potential /ixi in the TI defined with respect to the 
Dirac point. The last term in Eq. Q describes the tun¬ 
neling between points r of the SC and points x on the 
edge of the TI. The interface between the SC and the TI 
is assumed to be rough, which means that the thinnest 
regions of the interface give the highest probability for 
electrons to tunnel between the SC and the TI. We model 
these thinnest regions located at points xi by point con¬ 
tacts. The tunnel contact between the SC and the TI 
is doped with scalar and magnetic impurities with ran¬ 
domly oriented spins ^ = (5^^., Here ^ is 

the operator of the fc-th localized spin close to the point 
contact Xf, on the TI sample. The tunneling occurs via 
the virtual states of the scalar and magnetic impurities. 
The tunneling matrix amplitude is given by 


Ti(r,a;) = (5(r - a; ex) 

X 6(x — xg) 
e,k 


( 2 ) 


h+ 


Here, ex is a unit vector pointing along the TI edge writ¬ 
ten in terms of the coordinates of the SC. The normal 
(spin-flip) tunneling is parametrized by a complex am¬ 
plitude tk (Wfc) with scalar impurities contributing to the 
amplitude tk only. This model implies that there can be 
more than one magnetic or scalar impurity at the vicinity 
of the point contact. 

In the second model we consider the coupling of a bulk 
s-wave SC to a 2D TI via an array of single-level QDs, 
see Fig. 1(b). The Hamiltonian of the system is given by 


H 2 = Hbcs + Hti -I- Hd (3) 

+ ^ D] ■ T^^{Xj) + H.c. . 
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Here, Hd = (1/2) d] ■ + U, n,,;) + H.c. 

is the Hamiltonian of an array of single-level QDs at en¬ 
ergies €j > 0 and with amplitudes Uj of the Coulomb 
interaction on the QDs and Dj = 

is the electron Nambu operator on the jth QD. The oc¬ 
cupation number operators for spin-up and spin-down 
electrons on the jth QD are given by rij^^ = 

and rij^i = Tunneling occurs at points Vj and 

Xj on the SC and the TI, respectively, and is described 
by tunneling amplitudes tjp and tj^ 2 - 

Proximity-induced JttJs. We first discuss the model 
shown in Fig. Qa) and described by Eqs. Q and ([2]). 
We neglect the inverse proximity effect due to magnetic 
impurities. By integrating out the degrees of freedom 
of the SC and including contributions up to second or¬ 
der in the tunneling amplitudes we see that the equation 


(a) 


TI MI 
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FIG. 1: (Color online) Setups to generate a proximity-induced 
Josephson vr-junction in topological insulators (TIs). a) An 
s-wave SC (red) couples to a TI (grey) via an insulator doped 
with magnetic and scalar impurities (MI, magnetic insulator, 
blue). If the spin-flip tunneling rates are larger than the nor¬ 
mal tunneling rates superconducting gaps with opposite sign 
are induced in the TI samples, b) Instead of the MI the SC 
is coupled to the TI via an array of single-level QDs in the 
Coulomb blockade regime. The array of QDs is occupied with 
randomly oriented electron spins. 


of motion for the Green’s function g{x,x') of the TI in 
frequency space is given by 

g~^{x)-g{x,x') = S{x-x')-\- J dxiT.{x,xi)-g{xi,x') (4) 

with g~^{x) = iuj ittypck^dx and io the fermionic Mat- 
subara frequency. In leading order, the electron self¬ 
energy in the TI is given by 


'P{x,x\) = J d^r d^r' T|(r,x) • G(r 


r')-Ti(r',xi).(5) 


Here, G(r — r') denotes the Green’s function of the 
bare clean homogeneous three-dimensional SC defined 
by G“^(r) • G(r — r') = i5(r — r') with G“^(r) = iuj d- 
(k?d^/2m + fjjr^ — AscCr^ry. At vanishing relative dis¬ 
tance a solution to this equation is given 


= 0) = ■ 7 , ■ , - H > 

+ AL 


( 6 ) 


with V = the normal-state density of states per spin 
and pf the Eermi momentum in the SC. We adopt several 
assumptions to simplify Eq. ([^. Eirst, the distribution of 
impurities is assumed to be almost continuous and hence 
sums over impurities at discrete positions are replaced 
by integrals over impurity densities. Second, terms that 
are linear in the Pauli matrices (t“ vanish after averaging 
over the random orientation of the spins Third, at 
some fixed xg tunneling contributions from points xg> for 
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i' ^ I can be neglected. The contribution of these terms 
to the effective Hamiltonian can be incorporated in the 
chemical potential [40) . These assumptions imply that 


da::il](a;, Xi) • g{xi,x') 


iuj{r TTs) 


-A,e(r - 


9 {x,x') 


with the scattering rates 
r = TTt^nol 

k 

Ts = iTvnsS{S + l)Y,\ul\^ /i. 

k,a 


(7) 

( 8 ) 
(9) 


Here, uq is the concentration of point contacts that allow 
for spin-conserved tunneling, while ns is the concentra¬ 
tion of point contacts that allow for spin-flip tunneling. 
We note that = S{S + l)6£iiSkk>, with (...) 

meaning the average over random spin directions and S 
being the magnitude of the impurity spin. In particular, 
the average vanishes for different impurity spins. This 
implies that in the expression for the scattering rate Tg 
terms oc with k ^ k' vanish as well. The effective 

order parameter in the TI for lo ^ Age, T, Tg is given by 


Ar,rs~r-rs- (10) 


Interestingly, if Tg > T the effective order parameter can 
become negative. Such a situation naturally emerges if 
the tunnel contact contains a large number of magnetic 
and scalar impurities. At a particular point X£ the elec¬ 
tron tunneling amplitude via some magnetic impurity k 
is tk + k for some scalar impurity k'. 

We assume that \tk \ « \tk'\, while generally their relative 
sign is random. Thus, for many impurities the normal 
tunneling contributions in Eq. (|^ destructively interfere, 
so that J2k^k ~ 0. As a result, Tg > T can be realized 
and Ar.Fg becomes negative. 

Next we discuss the model of an s-wave SC coupled to 
a 2D TI via an array of QDs, as depicted in Fig. [^b) 
and described by Eq. ([^. We will work in the Coulomb 
blockade regime. Thus, we assume singly occupied QDs 
with the electron spin on the QDs being randomly ori¬ 
ented. In the limit of small tunneling amplitudes that 
couple the SC and the TI to the QD we use a Schrieffer- 
Wolff transformation [41] to map the Hamiltonian H 2 as 
given in Eq. onto a Hamiltonian Hi of the form as 
given in Eq. (W with tk = 0. The physical interpretation 
is that due to the large Coulomb interactions on the QDs 
only spin-flip tunneling of electrons through the dots is 
allowed [5| . From the discussion of the first model we can 
conclude again that Tg > 0, while T « 0. 

Thus, we see that in both models we obtain a JttJ in 
the tunneling region, be., the proximity-induced effec¬ 
tive superconducting order parameter in the helical edge 


states of the TI assumes the opposite sign compared to 
the one of the SC. 

Kramers pairs of Majorana fermions. In a TI 
proximity-coupled to an s-wave SC magnetic perturba¬ 
tions can be used to induce MFs |12|- However, in realis¬ 
tic scenarios the use of magnetic fields should be avoided 
since it acts detrimental on superconductivity and it is 
indeed not a necessary ingredient: In the absence of it 
Kramers pairs of MFs can emerge in nanowire systems 
that are coupled to unconventional SCs |2HH33]- Also se¬ 
tups using conventional SCs in proximity to nanowires 
[33H^ and in 2D [38] and to 3D TIs [39] were proposed. 
In particular, it was predicted that Kramers pairs of MFs 
appear due to JttJs in nanowires |33H35] or in 3D TI 
films [33] ■ In this section we make use of the JttJ mod¬ 
els introduced above and propose two setups (labeled by 
N = 1,2) that host Kramers pairs of MFs based on two 
2D TIs. As a major advantage both setups are accessible 
by current experimental techniques in TIs [3[3] and in 
JttJs based on magnetic insulators [3] and QDs [S]. Also, 
to reveal the non-abelian statistics the pairs can easily 
be moved by tuning a tunnel barrier between the TIs. 

We consider two TIs labeled by n = 1, 2. In the first 
(second) setup, edge states are of opposite (same) helic- 
ity and the chemical potentials are tuned to be opposite 
(to be the same) with /ii = —p 2 (mi = M 2 )) as illustrated 
in Fig. [2>(Fig . [^). Both TIs are brought into proximity 
to an s-wave SC. In the the first TI, the tunnel contact 
is doped with scalar and magnetic impurities with ran¬ 
domly oriented spins or, equivalently, an array of QDs 
with randomly oriented spins is used. As shown above, a 
TT-junction emerges and the proximity gap in the first TI 
acquires the opposite sign to the bulk SC, — Ai < 0. The 
tunnel contact between the SC and the second TI does 
not contain a spin-flip source. Thus, the corresponding 
order parameter is of the same sign as in the SC, A 2 > 0. 
The induced superconductivity in the nth TI of the Ath 
setup is described by the Hamiltonian 




dx 


Ll Rl - U + H.c. 

( 11 ) 


in the basis 4>„(x) = (R„(a:), L„(a:), r 1 (x), L))(x)). For 
the first (second) TI of the first setup we have intro¬ 
duced slowly-varying spin-up (spin-down) right-mover 
Ri(x) [R 2 (x)] and spin-down (spin-up) left-mover Li(x) 
[L 2 (x)] fields defined around the Fermi points dzkp. In 
the second setup R 2 (x) [L 2 (x)] is the spin up (spin down) 
mode, see Fig. [^ 

The two TIs are coupled via a tunnel barrier placed in 
the region 0 ^ x ^ L as shown in Fig. [^a) Neglecting 
the fast-oscillating terms [45] , we find that the tunneling 
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FIG. 2: (Color online) (a) Setup hosting Kramers pairs of 
MFs. Two TIs (grey rectangles) are placed on top of an un¬ 
derlying s-wave SC (red) such that proximity superconduc¬ 
tivity is induced in both TIs. Importantly, the first TI is 
coupled through a magnetic insulator (MI, blue) resulting in 
the JttJ. The tunnel barrier (TB, yellow) between the edges 
of two TIs extends from x = Q to x = L. One Kramers pair 
of MFs 7 i ,2 [ 73 , 4 ] (purple) is localized at the a; = 0 [ 2 ; = L] 
end of the TB and can be manipulated by tuning the length 
L of the TB. The spectrum of two pairs of TI edge modes 
is considered in (b) for the first setup and in (c) for the sec¬ 
ond setup, (b) Edge modes of the same TI are coupled by 
proximity-induced pairing amplitudes Ai < 0 and A 2 > 0, 
resp. The chemical potentials are opposite for the two TIs, 
fii = —p 2 - The helicities of the edge states are opposite (in¬ 
dicated by the coloring in red and blue). The tunneling (t) 
couples a right-moving state in the hrst TI to a left-moving 
state in the second TI, and vice versa, (c) The two TIs have 
the same chemical potential, pi = 1 x 2 , and the same helicities. 
The tunneling (t) couples a right-[left-] moving state in the 
first TI to a right-[left-] moving state in the second TI. 


Hamiltonian in the first setup is given by 

= dx[e^^ (r^Li - LiR|) (12) 

+ e-^^ (l^Ri -RiLJ) -bH.c.], 
while in the second setup by 

= dx[e^^ (r^Ri - RiRJ) (13) 

-|- e -j- H.c.]. 

Here, t {(p) is the tunneling amplitude (phase) between 
two TIs. The choice of 4> ensures TRI. The total Hamil¬ 
tonian is + Hic,n) + where the 

kinetic part Hxi^n = Hpi is identical for both TIs and 
was introduced in Eq. Q. The tunneling phase can 
be removed from the total Hamiltonian by a suitable 
gauge transformation [33]. In both setups, we find that 
a topological phase transition accompanied by the bulk 
gap closing and reopening occurs at the point 

(14) 


In the second setup, there is an additional constraint 
Ai 7 ^ A 2 . If t > v^AiA 2 , there is one Kramers pair 
of MFs localized at the interfaces at cc = 0 and one at 
X = L |33| . The localization lengths are inversely propor¬ 
tional to the gaps opened at the Fermi points |33|- Thus, 
in regions with no tunnel coupling between the TIs the 
localization lengths are the superconducting coherence 
lengths = ft-vp/An, while in regions with t > A 1 A 2 
they are given by 


= 2hvF/W(Ai-A2Y + Af^- Ai - A 2 ) 



2hvF 

A2r± aVT^l + A2)2 - 4t2 ■ 


(15) 

(16) 


Superscript (1) [(2)] corresponds to first [second] setup 
and 5ft means the real part of a complex number. In both 
setups we assume that the length of the tunnel barrier 
L is much longer than the localization lengths ^max = 
max{^i,^ 2 ,^^^^} (Cmix = max{^ 2 ,'?i^^} for Ai > A 2 and 
Cmax = max{^i, for Ai < A 2 ). Hence the wavefunc- 
tions of the MFs localized at the two different interfaces 
do not overlap and can be considered as independent. If 
L is comparable or shorter than the localization length of 
the MFs they hybridize into two complex fermionic state 
whose energies are non-zero in general [451146] . Tuning L 
by underlying gates allows one to arbitrarily control the 
position of the MFs along the TI edges. To give a numer¬ 
ical estimate of the localization length we assume that 
the induced gaps are given by Ai = 0.1 meV, A 2 = 0.2 
meV, and the tunnel coupling is set to t = 0.2 meV. 
In a InAs/GaSb (HgTe/CdTe) TI the Fermi velocity is 
given by vp = 4.6 x 10"* ms“^ [5] (5.5 x 10® ms“^ [ 11 ]). 
This gives a localization length of the order of 0.5 /xm 
(5 /xm). Finally, we emphasize here that the opposite 
relative sign in front of the proximity induced gaps of 
the two edges in Eq.(ll) is an important ingredient to 
generate the Kramers pair of MFs. If this relative sign is 
the same there exists no topological phase. For illustra¬ 
tive phase diagrams, see [44] . 

Conclusions. We have proposed and studied two se¬ 
tups to realize a proximity-induced JttJ in a TI in the 
presence of TRI. Both setups rely on the tunnel cou¬ 
pling of a TI sample to an s-wave bulk SC either via a 
layer of scalar and magnetic impurities with randomly 
oriented spins or via an array of QDs each of which is 
occupied by a randomly oriented spin. We have seen 
that if in either case spin-flip tunneling dominates over 
normal tunneling a xr-junction emerges. The randomly 
oriented spins ensure that there is effectively no break¬ 
ing of TRI. As an application we have demonstrated how 
such proximity-induced tt-J unctions can be used to gener¬ 
ate and manipulate Kramers pairs of MFs in edge states 
of tunnel-coupled TIs. 
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In the Supplemental Material, we derive the MF wave- 
functions for two models introduced in the main text. 


Energy spectrum 


if t > Ai = A 2 . We now assume that Ai ^ A 2 and con¬ 
firm that Eq. (S3) defines a topological phase transition. 
This means that there should be MFs localized at the 
boundary between two space regions with t > A 1 A 2 
and t < yfAiA 2 . 


We find that the bulk spectrum of the Hamiltonian 
from the main text is given by 


Wavefunctions and localization lengths of MFs 


E- 


l,s,± 


(k) =. 


{hvFkY + (^A+ ± ^A^_+ 


-11/2 


(SI) 


where k is the momentum in the TI, and s = ±1. Sim¬ 
ilarly the bulk spectrum of the Hamiltonian from 
the main text is given by 


E 2 ,s,±{k) = s 


{hvpkf P A\ + A'i+f ± 2^/W{k) 

(S2) 


with W{k) = {hvpkYt^ + A^(A^ -|- 1 ^). Here, we also 
introduced the notations A± = (Ai ± A 2 )/ 2 . Both spec¬ 
tra Ei^s,±{k) and £ 2 , 3 ,±{k) are twofold degenerate as ex¬ 
pected for time-reversal invariant systems. 

We find that the spectrum is gapless at fc = 0 if 


t = v'AiA 2 , (S3) 


and is gapped otherwise. Here, for setup N = 2 we need 
the additional condition that Ai 7 ^ A 2 . Also the spectral 
gap for the setup N = 2 closes at some finite momentum 


The operator defining a MF, which is a zero-energy 
bound state, is generally given by = (yj^^)^ = 

J2n=i 2 I V'i^^(^) ■ ^n{x) with the wavefunction (vec¬ 
tor) ’ 


{if, 


{N)\T 
n-d ) 


[x) 


f^Z^{x) 

{&*{x) 

yialVyix)/ 


(S4) 


for some complex-valued functions fZj {x) and gZ\x). 
The index j = 1,2 distinguishes between two MFs be¬ 
longing to the same Kramers pair. The form of these 
functions is different for different setups. 

Without loss of generality, we focus below on the left 
interfaces at which the tunneling amplitude jumps from 
t = 0 ata;<0tot = to> VA 1 A 2 for a; > 0. 

First setup. We find that for the first setup the inter¬ 
face hosts a Kramers pair of MFs given by 


Z ^ + Sn 2 e e if x > 0 

gT if x < 0 

fZ = (-l)”(5ia)*> 9nl = (-l)”-'(/S)* 


-ifnl = {gnZ* = 


(S5) 


with the localization lengths given by 

^(1) = hvF/{\J A'^_+tl- A+) 

~ F!^n- 

In Fig. Fla) the localization length is plotted 
for different valus of to in color scale versus Ai and 


A 2 . Note that the solutions for given N are orthogo¬ 
nal, ip^i ■ {ifZz)'^ — i'hns forming a Kramers pair. 

The localization length of the MF is ^en by ^max = 
max{^(^i,^ 1 ,^ 2 } and is plotted in Fig. Fl^)- 

Second setup. The interface at x = 0 of the second 
setup also hosts a Kramers pair of MFs. For A_ > 0 the 
Kramers pair of MFs is given by 
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f(2) 


1(5. 


n2 


f(2) _ „ikpx 

Jn,l — ^ 


„—i 


2y/A^^-t 


- +iSu 2 il + e 

_ / ^ to_sin(fc(^ 4 _ 1 _ ,-A „ 


e< 2 ) 


^±Z^SHlp-“/«+’ - ^ + + V^+-*V-x/;<^^A _ S^ije "^4 )toe-** 


,( 2 ) 


e 2 


—55„2 e* 2 e “/^2 


sgl=.(/S)-. /S = -toar 


„( 2 ) _ I/•(!>(. 

i/n,2 ~ Wn,lJ 


2 ^AF«^ 


cos{k^‘^^x) ] e *2 ] e 


( 2 ) 


if a; > 0, A+ > fo 
if a; > 0, A+ = to 
if a; > 0, A+ < to 
if a; < 0 

(S7) 


while for A_ < 0 it is given by 


f(2) _ „ikpa 

Jn,l — ^ 


iSn 




( 2 ) 

+ —^ 


-x/i 


( 2 ) 


toe* 2 


2y/A^-t 


— <5nl 


A±-yAF7| /^(2) _ A+ + yA^ /5^2) 


2 .yAF* 


5^2) 




Ki ( + ^os(k»>x)^ 'F-™r 


.<5nl 


e -*2 e“/«i 


( 2 ) _ „ . ,( 2 )x* ,( 2 ) _ ( 2 ) ( 2 ) _ . ,( 2 ).* 

UnA — ‘Urt.l/ J an,2 ~ KUn,!) ) ^71,2 ~ Un.ll 


e “*2 if a: > 0, A+ > to 
if a: > 0, A+ = to 
if a; > 0, A-). < to 
a X <0 

(S8) 


Again, the two wavefunctions for given N are orthogonal 
and thus forming a Kramers pair. We have defined the 
wavenumber 


= ^J\Al-tl\/hvF 

and the localization lengths 

p^{2) _ _ _ flVF __ 

|A_|±5i^A2^-t2 


(S9) 


(SIO) 


(Cmax = max{^i,and is plotted in Fig. ^e). The 
fast-oscillating factors have explicitly been restored in the 
wavefunctions. We see that for Ai = A 2 the MF wave- 
function is delocalized. In the limit when L ^ ^inlx the 
interfaces at a; = 0 and x = L can be considered as in¬ 
dependent and a calculation of the MF wavefunctions at 
X = L can be performed analogously. 

Finally for setup N the tunneling phase (f) can be ab¬ 
sorbed into a redefinition of the electron operators by the 
gauge transformations 


( 2 ) 

For different values of to the localization lengths F are 
displayed in Fig. % c) and (d) in a color scale plot as 
a function of Ai and A 2 . For to < A+ the spectral gap 
closes at zero momentum while for to > A+ it closes at 
some finite momentum. For A_ > 0 (A_ < 0) the local¬ 
ization length of the MF is given by ^max = max{^ 2 , 


R„ 1 -^ exp 

L„ 1 -^ exp L, 


(Sll) 




































Figure SI: (Color online) (a) Phase diagrams and color scale plots of the localization length for a: > 0 of the Kramers 
pair of MFs in the first setup as a function of the superconducting gap parameters Ai ,2 and the tunneling amplitude to- Here, 
increases from blue, through yellow, to red; vp = 4.6 x 10"^ ms“^ in an InAs/GaSb TI [5]. The curve to = \/AiA 2 
(dashed) seperates the topological phase (T, colored) and the non-topological phase (NT, uncolored). At the phase boundary 
the localization length is divergent. For a: < 0 the localization lengths are given by the superconducting coherence lengths 
Cl ,2 = hvF/^ 1 , 2 - (b) Same as in (a) but for Cmix = max{Ci,Ca,The curves Ci ,2 = a-nd Ci = Cz (solid) seperate 
regions where Cmlx is given respectively by CijCz or (c) Same as in (a) but for the localization length for a: > 0 in 
the second setup. Along the line Ai = Aa (dashed) the localization lengths are divergent, (d) Same as in (c) but for the 
localization length (e) Same as in (d) but for Cmlx = max{C2,C-^} T Ai > A 2 and Cmlx = max{Ci,C-^} T < A 2 . The 
curves Ci ,2 = C- (solid) seperate regions where Cmax is given respectively by C 11 C 2 or . 































































